Statistical model of dephasing in mesoscopic devices introduced in the scattering 

matrix formalism 
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We propose a phenomenological model of dephasing in mesoscopic transport, based on the intro- 
duction of random phase fluctuations in the computation of the scattering matrix of the system. A 
Monte Carlo averaging procedure allows us to extract electrical and microscopic device properties. 
We show that, in this picture, scattering matrix properties enforced by current conservation and 
time reversal invariance still hold. In order to assess the validity of the proposed approach, we 
present simulations of conductance and magnetoconductance of Aharonov-Bohm rings that repro- 
duce the behavior observed in experiments, in particular as far as aspects related to decoherence 
are concerned. 



i 

-a 
c 

o 
o 



> 
00 

m 
O 
-i— > 



i 

C 

o 
o 



I. INTRODUCTION 

Phase coherence of the electron wave function has 
a fundamental influence on the transport properties of 
mesoscopic devices^ and is at the basis of several phenom- 
ena ranging from interference effects (such as Aharonov- 
Bohm oscillations) to weak localization (WL) 2 ^ and uni- 
versal conductance fluctuations (UCF)4£ 

The fundamental quantity used to express the degree 
of coherence in a system is the so-called phase coherence 
length (or "dephasing" length) that is typically esti- 
mated on the basis of WL experiments in semiconductor 
heterostructures^ Si MOSFETs 8 , metal conductors^ or 
of interference experiments in devices such as Aharonov- 
Bohm (AB) rmgsiGiii 

Mesoscopic physics deals with devices whose size is 
smaller or comparable to 1$ and therefore often operate 
in an intermediate regime between coherent transport, in 
which the phase information is fully preserved, and inco- 
herent transport. The main phase-breaking mechanisms 
are due to interaction of electrons with other electrons, 
photons, phonons, and defects such as magnetic impu- 
rities or to other kinds of phase-randomizing interaction 
with the environment. 12,13 

Therefore, it would be very useful to have a unique 
formalism capable to include an arbitrary degree of de- 
phasing in the evaluation of the transport properties of 
a system, and to allow a seamless transition between the 
coherent and the fully incoherent limits. 

In the case of interfering paths, an "ad-hoc" random 
term can be added analytically to the difference between 
the phases accumulated in the two paths. When generic 
devices with two or more leads are considered, two main 
phenomenological models are available for including a 
partial degree of dephasing in the transport model: i) 
insertion in the device of an additional "virtual" voltage 
probed that can be also taken into account by properly 
adjusting the two-terminal conductanc e) 15 ! 16 ; or ii) ad- 
dition of an imaginary part to the Hamiltonian in the 
device region piLiLiS I n case i), the seeming drawback 
of spatially localized decoherence can be overcome either 



by introducing an adequate number of virtual probes in 
different points of the device region^ or by considering 
the limit of a voltage lead that supports an infinite num- 
ber of modest 5 - In case ii), the carriers absorbed by the 
imaginary term have to be re-injected into the conductor 
in order to ensure current conservation. 

An additional method to treat dephasing consists 
in including a stochastic absorption in the scattering 
descriptionSLSSiS through the insertion of an attenuation 
factor in the free propagation region. Also in this case, 
continuity of the probability density current requires that 
absorbed electrons are re-injected. 

Dephasing due to the environment can be modeled by 
two equivalent approaches^ One focuses on the changes 
that the wave function induces on the state of the envi- 
ronment, and was adopted, for instance, to simulate elec- 
tron conduction interacting with dynamic impurities. 24 
The other addresses the phase accumulated by the inter- 
fering waves as a statistical process. In this paper, we 
adopt the latter perspective and propose a phenomeno- 
logical model of decoherence, that treats dephasing as 
a distributed phenomenon in the device region, ensures 
the conservation of current density and allows to evaluate 
the local density of states. We consider the stochastic be- 
havior of the dephasing process and adopt a Monte Carlo 
averaging procedure to extract the electrical and micro- 
scopic properties of the system. We are able to vary 1$ 
and gradually move from a coherent to a totally incoher- 
ent transport regime. The model is described in Sec. ITT1 
and is applied in Sec. IHII to evaluate the decoherence on 
the conductance and magnetoconductance of an AB ring. 



II. DEPHASING MODEL 

We include our model for dephasing in the scattering 
matrix formalism for the computation of the device con- 
ductance G. The conductance of a generic structure is 
related to the transmission probability matrix T = t't by 
the Landauer-Biittiker formula G = ge 2 /h J2 n m T nm r^ 
where t is the transmission matrix, g is the spin degen- 
eracy factor, e is the elementary charge, h is Planck's 



2 



constant and n, m run over all transverse modes con- 
tributing to transport. 

The transmission matrix is obtained by computing the 
scattering matrix (S-matrix) of the device^ If we subdi- 
vide the domain along the transport direction x in several 
slices, one for each grid-point j = 1, . . . , N x in the x di- 
rection, the wave function of the electron in the j-th slice 
(xj < x < Xj + i) can be written as 

^( X , y ) = ^^^(a j , n e ik ^+b j , n e- ik ^) , (1) 
„ V\kj,n\ 

where Xj,n{y) is the n-th transverse eigenvector of the j- 
th slice with eigenenergies Ej >n and the longitudinal wave 
vector kj >n is related to the total energy E by the con- 
dition E = Ej. n + h 2 k 2 n /2rrij. The coefficients aj.„ and 
bj. n are obtained by imposing the continuity of the wave 
function and of the probability current density at the 
interface between the j-th and the (j + l)-th slice. The 
scattering matrix Sj links the incoming and the outgoing 
coefficients: 




where cij (bj) is the column vector of all aj .„'s (& Jin 's) 
for n = 1, . . . , A^modc, and N mo d c is the total num- 
ber of modes considered in the system. The compo- 
sition between two adjacent scattering matrices^, Sj 
and Sj+i, gives the matrix Sj ® Sj+i which links aj, 
bj, Oj+2, and 6j+2- In order to compute the scatter- 
ing matrix of the complete device we have to compose 
the matrices of all slices according to well known rulesi^l 
S = Si <g> •• • <g> Sj ® •• • ® S Nx -i- 

We model the effect of decoherence as a random vari- 
ation of the phase accumulated by each mode in each 
slice into which the device has been divided. In the ab- 
sence of dephasing, mode n accumulates in slice j a phase 
kj,n(%j+i ~ Xj); in the presence of dephasing it accumu- 
lates a phase kj^ n (xj+i — Xj) + A<pj iT1 , where ^4>j, n is 
a random term obeying a Gaussian probability distribu- 
tion with zero average and standard deviation aj, that 
depends on the thickness of the slice (a^+i — Xj) and on 
1$, as we shall show. 

For a random choice of all A</>j jn 's, for j = 1, . . . , N x — 
1, n = 1, . . . , N mo de, we can compute an "occurrence" 
S of the scattering matrix of the system. We take into 
account the probabilistic nature of dephasing, and there- 
fore transport properties are obtained following a Monte 
Carlo averaging procedure over a sufficiently large en- 
semble of random occurrences. For typical devices, in 
order to obtain stable and "smooth" averages, we need 
to consider an ensemble of about a hundred occurrences. 

For the purpose of clarity, we have described the case of 
a two-terminal devices. However, the method can be ap- 
plied without any variation for the computation of many- 
terminal scattering matrices. 

Using some algebra it is straightforward to verify the 
unitarity of any S: adding the random phase term to the 



scattering matrix of the j-th slice corresponds to sub- 
stitute the coherent matrix Sj with Sj = Sj ® S™ ndora , 
where Sj andora is a scattering matrix in which the reflec- 
tion matrices are zero, and the transmission matrices are 
diagonal, their n-th element (n = 1, . . . , N mo d e ) being 
exp(iA(j) hn ). 

It is easy to verify that 5'j andom is unitary by construc- 
tion; since composition of unitary scattering matrices 
provide a unitary scattering matrix, each S is unitary. 

The physical reason for unitarity of S is the conserva- 
tion of the incoming current, whereas the time reversal 
symmetry in the presence of a magnetic flux $ implies 
the validity of the Onsager-Casimir relations^ 

T M ($) = T w (-$) , R PP (<i>) = Rppi-*) , (3) 

where the labels q and p denote the leads of the sys- 
tem and T pq is the total transmission probability from 
lead p to lead q (summed over all modes), and R pp is 
the total reflection probability at lead p~ Once again, 
each Sj andoin is symmetric and independent of the mag- 
netic field and therefore obeys 0; it is now sufficient to 
observe that a composition of matrices obeying 10) still 
provides a matrix that obeys Onsager-Casimir relations. 

Let us consider a traveling plane wave that loses phase 
coherence as it propagates, but conserves its modulus. 
One possible description of such situation is to write the 
wave function as the sum of a coherent component whose 
amplitude decays exponentially with propagation for a 
length I as e ^ V 2 ^ an d f an incoherent component to- 
tally uncorrelated with the former, that ensures conser- 
vation of the wave function modulus. Another possible 
description is to add to the phase of the traveling wave 
function after a length I a random term, with Gaussian 
distribution, zero average and standard deviation a. 

In order to derive the relation between and a, we 
consider the case of wave interference. First, let us con- 
sider two coherent wave functions ipi and ^2 of amplitude 
unity, obtained for example with a beam splitter. We let 
them interfere again after both propagate along paths of 
length I. In terms of the former description the amplitude 
of the interfering pattern is 

+ HLx - l^i + iMmin = 4expH/W. (4) 

On the other hand, if we write the same two wave func- 
tions with the latter description, they have amplitude 
unity and phases containing additional random terms <^ 
and <ftft, respectively, that are uncorrelated, and obey a 
Gaussian distribution with average zero and standard de- 
viation a. The amplitude of the interfering pattern, in 
this case, is 

(| e ^R + e^R| 2 )-(| e ^a_ e #R|2) = 4(cos(^-^)) = 

(5) 

where angle brackets mean statistical averaging, and the 
last equation has been obtained using the fact that </>r — 
§\ is a Gaussian variable of average zero and variance 
2a 2 , and (cos0) = J d<\> coscj) exp[— (j) 2 /4a 2 } /Vina 2 = 
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By comparing @ and we obtain a 2 = l/l<p, that, 
if we consider each single slice in which the structure is 
partitioned, means 




III. SIMULATIONS 

In this section we use our proposed model for de- 
phasing to investigate the effect of decoherence on UCF 
and Aharonov-Bohm oscillations in mesoscopic rings, for 
which experimental results are available in the literature. 

The Aharonov-Bohm ring (shown the inset of Fig.^l is 
broadly used to perform phase coherence measurements 
because it provides the possibility to obtain WL, UCF, 
as well as pure interference effectsi As far as this as- 
pect is concerned, analytical predictions of the dephas- 
ing rate are available^ that agree very well with the 
experiments^ 

UCF appears when an external parameter that alters 
the potential profile of the structure is varied. Indeed, 
such conductance fluctuations are obtained in experi- 
ments by varying the Fermi level Ep of the electrons 
through the voltage on the back gate or on a top gate. 
The typical amplitude of conductance fluctuations does 
not depend on the sample size or on the degree of disor- 
der and is of the order of the conductance quantum 2e 2 /h 
in a purely coherent transport regime4^*i& If this is not 
the case, decoherence smears out fluctuations restoring 
a staircase when G is plotted versus Ep or versus the 
gate vo ltag e, as shown for example in the experiments of 

Ref . lillliT 

We have simulated a symmetric AB ring structure 
defined by etching on a GaAs/AlGaAs heterostructure. 
With reference to the inset of Fig. ^ the internal radius 
of the ring is 350 nm, the external radius is 630 nm, the 
width of the leads is 200 nm. Bright regions correspond 
to a potential energy of eV, dark regions to 0.2 eV. 

The thin line in Fig. ^ is the computed conductance as 
a function of Ep for fully coherent transport, while the 
thick line is the conductance corresponding to a dephas- 
ing length 1$ = 0.3 /im. Results are obtained by averag- 
ing on 100 random occurrences. The fluctuations clearly 
present in the coherent regime are evidently smoothed 
out as decoherence is introduced. 

In addition, as a consequence of the introduction of 
decoherence we are able to observe the loss of the univer- 
sality in the conductance behavior, that now presents a 
series of non-integer plateaus. Such behavior is typically 
observed in experimentsiS^S^i as is due to backscatter- 
ing phenomena and to scattering at the interface between 
ID and 2D electron gas. However, some degree of deco- 
herence, which is always present in experiments, is re- 
quired to clearly reproduce the phenomenon with simu- 
lations. 
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FIG. 1: (Color online). Conductance of the AB ring as a 
function of the Fermi level of electrons: Completely coher- 
ent regime (thin solid line), partially incoherent regime cor- 
responding to l,), = 0.3/im (thick solid line). Results are av- 
eraged on 100 random occurrences. The arrows indicate the 
energies at which a new conducting channel in the lead is 
opened. Inset: the AB ring potential used in our simulations. 
The internal radius R a is 350 nm, whereas the external radius 
Tib is 630 nm. The width of the branches W that connect the 
ring to drain and source is 200 nm. Go is the conductance 
quantum 2e 2 /h. 

We want to emphasize that dephasing introduced by 
our model has a very different effect on the device conduc- 
tance than energy averaging due to a finite temperature 
of the system. In Fig. we plot the conductance of the 
same ring in Fig. ^ versus Ep and compare it with the 
thermal-averaged conductance computed in the case of 
fully coherent transport: 

G{T) = -jdEG{E) d J^l, 

where f(E,T) is the Fermi-Dirac occupation factor. As 
can be seen, the coherent conductance at different tem- 
peratures does not exhibit the non-integer plateaus previ- 
ously observed and has a behavior qualitatively different 
from partially coherent conductance. 

Another way to verify how dephasing influences the 
transport properties of the ring structure is represented 
by the study of magnetoconductance. In our code we 
have added the effects of an external magnetic field 
B = (0, 0, B) perpendicular to the propagation plane 
xy. We adopt the transverse gauge A = (—By, 0, 0) for 
the vector potential A as described in Ref. |33- Due to 
the AB ring geometry the phase difference of wave func- 
tions propagating along the two branches depends on the 
magnetic field as J(p — eA)/h ■ dr, generating the mag- 
netoconductance oscillations. The oscillation period can 
be equal to the quantum flux h/e or to the submultiplcs 
h/ne when coherent backscattering is present and the 
electron turns around the ring more times. As expected, 
decoherence suppresses the amplitude of magnetoconduc- 
tance. Results for the AB ring geometry of the inset of 
Fig. n are shown in Fig. [3 where it is possible to appre- 
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FIG. 2: (Color online). Comparison between the conductance 
of the AB ring at K computed with 1^ — 0.3 /im (thick 
line) and fully coherent conductances computed at finite tem- 
perature (thin lines). For clarity of presentation each line is 
shifted by one unit of conductance 



ciate the transition from a coherent transport regime to 
an only partially coherent one as Ia, is decreased. 
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FIG. 3: (Color online). The conductance oscillations of the 
AB ring shown in the inset of Fig. Q as a function of the 
magnetic field, for dephasing length /J, = oo, 1^ — 5 fim, = 
1 /im, — 0.3 fim. G is averaged over 100 simulation runs. 
For clarity of presentation, each line is shifted by one conduc- 
tance unit. 

In Fig. ^ we show a comparison with experimental re- 
sults presented by Hansen et alm Q for a symmetric AB 
ring with internal radius 280 nm, external radius 560 nm, 
and wire width 100 nm. On the left we show the Fast 
Fourier Transform (FFT) of experimental magnetocon- 
ductance oscillations (Fig. 2 of Ref. IT(ih for different val- 
ues of temperature, while on the right we show the same 
quantity computed with our model for different values of 
the dephasing length. In both cases all frequency com- 
ponents are damped by decoherence. At small values of 

only the first peak corresponding to the h/e frequency 
is clearly visible. 

It is important to verify whether the dephasing length 
microscopically introduced by our model through Eq. © 
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FIG. 4: (Color online). Left: Fast Fourier Transform (FFT) 
of the AB experimental oscillations measured by Hansen et al. 
(Fig. 2 of Ref. lToF) at different temperatures. Right: the same 
FFT obtained with our simulations for values of the dephasing 



length l\ = oo, \\ 



fim. 
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1 /mi, Zj, = 0.3 /im. In the 



simulation results, the DC component has been removed. 



agrees with the value that can be extracted from the 
electrical properties of the device. Following Ref. we 
assume that the amplitude A n of the h/ne oscillation can 
be written as 



A n cx e ~ nL ^ 



(7) 



where L is the circumference of the ring. Such assump- 
tion is confirmed by the experimental results shown in 
Fig. EJi (Fig. 2 of Ref. [F 
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FIG. 5: Left: Measured oscillation amplitude for the h/e 
and /i/4e frequencies (Fig. 2 of Ref. HfTl . Right: Amplitude of 
the h/e and h/3e oscillations plotted as a function of the de- 
phasing length Ij,. Each dot corresponds to a different Fermi 
energy of the propagating electrons. 

The FFT of the simulated oscillation amplitudes for 
different n exhibit an exponential dependence on 1$ as 
shown in the semilog plot of Fig. 03 for the cases n—1 
and ro=3. Both the slopes of the h/3e and the h/e os- 
cillations, according to Eq. (JJJ, are consistent with the 
nominal value provided by Eq. 10. Results for n > 3 
are not reliable in the whole range of 1$ due to numerical 
fluctuations and therefore are not shown. 
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FIG. 6: (Color online). Partial density of states for different 
dephasing lengths in the AB ring of Fig. 1 when only the first 
mode is populated. The magnetic field B is 55 mT and 85 mT, 
corresponding to a maximum and a minimum of conductance 
(see Fig.EJ. 



Finally, we emphasize the possibility of a microscopic 
description of the effects of decoherence in the system. In 
Fig. El we show the partial local density of states corre- 
sponding to one mode injected from the left p(x, y, E) oc 
\tp(x, y, E) | 2 for three different values of l<f, and two differ- 
ent values of the magnetic field B, corresponding to the 
cases of maximum constructive interference (B = 55 mT) 
and maximum destructive interference (B = 85 mT), as 
can be seen in Fig. 3. Also in this case the density of 
states is obtained by averaging over 100 Monte Carlo 
runs. 

For 1^ — 5 /im transport is almost fully coherent, a 
clear pattern of nodes forms in both branches, and in 
the output lead we have maximum modulation of the 



density of states as a function of B. For smaller dephas- 
ing lengths the stationary wave pattern in the branches 
smooths out. In particular, for 1$ = 0.3 /im, when the 
interference pattern is almost destroyed, as can be seen 
in Fig. 3, the density of states is quasi-constant in the 
branches and in the leads. 

IV. CONCLUSIONS 

In this paper we have presented a phcnomcnological 
microscopic model for the simulation of dephasing in 
mesoscopic devices. The stochastic nature of the dephas- 
ing process is taken into account with a Monte Carlo 
methodology used to extract average conductance, mag- 
netoconductance and density of states. We have shown 
that the proposed method ensures the physical validity 
of each occurrence of the scattering matrix. 

Here we want to underline the fact that the proposed 
method provides a unique description applicable to sys- 
tems with an arbitrary degree of dephasing. This repre- 
sents the main advantage of the proposed method, since 
common methods for determining transport properties 
of generic devices consider only the limit of completely 
coherent transport (with scattering matrix or recursive 
Green's functions techniques) or the limit of fully inco- 
herent transport (with semiclassical approaches). 

We have also shown that such method allows to recover 
experimental results observed in Aharonov-Bohm rings, 
when a certain degree of decoherence is always present 
and responsible for some typical features, such as the 
non-integer conductance plateaus at zero magnetic field. 

We believe that the proposed model can be very useful 
in understanding the effect of dephasing on the trans- 
port properties of mesoscopic devices, and enables to ac- 
curately reproduce experimental results with numerical 
simulations. It can also have a significant effect in as- 
sessing the effect of dephasing on the noise properties of 
nanoscale devices. 
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